Entanglement between distant qubits in cyclic XX chains 
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We evaluate the exact concurrence between any two spins in a cyclic XX chain of n spins placed 
in a uniform transverse magnetic field, both at zero and finite temperature, by means of the Jordan- 
Wigner transformation plus a number parity projected statistics. It is shown that while at T = 
there is always entanglement between any two spins in a narrow field interval before the transition to 
the aligned state, at low but non-zero temperatures the entanglement remains non-zero for arbitrarily 
high fields, for any pair separation L, although its magnitude decreases exponentially with the field. 
It is also demonstrated that all associated limit temperatures approach a constant non-zero value in 
this limit, which decreases as L~ 2 for L <n but exhibit special finite size effects for distant qubits 
(L ~ n/2). Related aspects such as the different behavior of even and odd antiferromagnetic chains, 
the existence of n ground state transitions and the thermodynamic limit n — > oo are also discussed. 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 



I. INTRODUCTION 

Quantum entanglement denotes those correlations 
with no classical analogue that can be exhibited by com- 
posite quantum systems and that constitute one of the 
most fundamental features of quantum mechanics. It is 
considered an essential resource in the field of quantum 
information [l[ , where it plays a key role in various quan- 
tum information processing tasks such as quantum tele- 
portation Q and quantum cryptography |3[. It is also 
playing an increasingly important role in condensed mat- 
ter physics, providing a new perspective for understand- 
ing quantum phase transitions and collective phenomena 
in strongly correlated systems [J, H, IE • 

In particular, there has been considerable interest in 
investigating entanglement in quantum spin chains with 
Heisenberg interactions @, Q , since they provide a scal- 
able qubit representation apt for quantum processing 
tasks [Hi EH whi ch 

can be realized in diverse physi- 
cal systems. Studies of thepairwise entanglement in 
the Ising and XY models [J . [a 12| and in the isotropic 
Heisenberg model [H, EH, EH, El | at zero and finite tem- 
perature and in a transverse uniform field, as well as in 
diverse XX, XY and XYZ models for two or a small 
number of qubits [13, El, El, [H, have been made. 
An important result is that the entanglement range may 
remain finite at a quantum phase transition, limited for 
instance to first and second neighbors in the Ising model 
0, 0], in contrast with the behavior of the correlation 
length, which diverges at these points. Global thermal 
entanglement has also been studied [22], showing that 
limit temperatures for pairwise entanglement arc lower 
bounds to those limiting entanglement between global 
partitions. A fundamental result for finite systems is 
that there is always a finite limit temperature for en- 
tanglement, since any mixed state becomes completely 
separable if it is sufficiently close to the full random state 

mm. 

In this work we analyze the entanglement between any 
two spins in a cyclic chain with nearest neighbor XX 
coupling in a transverse magnetic field (control param- 



eter) by means of an exact analytic treatment valid for 
any spin number n and pair separation L, based on the 
Jordan- Wigner mapping and the use of number parity 
projected statistics for T > 0. Recent studies in XX 
chains have focused either chains with a small number of 
spins [It], [HI HH , where results were obtained through 
direct diagonalization, or open chains at zero tempera- 
ture and field [l5|. We will show that the XX model 
offers very interesting properties such as entanglement 
between any pair (full range) in a finite field interval just 
before the critical point at T — 0, which subsists for large 
fields at low but non-zero temperatures T < Tr,. More- 
over, limit temperatures Tj, approach a non-zero limit 
for large fields, for all separations L. It also displays n 
ground state transitions at analytic field values, entailing 
a stepwise variation of the entanglement range suitable 
for its use as an entanglement switch. Let us mention 
that XX chains have also been employed for entangle- 
ment teleportation [2a |. 

Section II describes the formalism for evaluating the 
exact concurrence between arbitrary sites both at zero 
and finite temperature. Section III describes the main 
physical results, including the ground state transitions 
and concurrence both in ferro- and antiferromagnetic sys- 
tems, and a detail study of the limit temperatures for 
entanglement. Conclusions are drawn in IV. 



II. FORMALISM 

We consider a cyclic chain of n spins with nearest 
neighbor XX coupling. The Hamiltonian reads 



h = bs*-v-£(s*s* +1 + s y jS « +1 ) 



(la) 



= b s z - 5«E(4 s 7 + i + s 7+i4) . ( lb ) 

where s*' v ' z are the spin components (in units of H) at site 
j, = sj±is|, S z — 5j? =1 s z - is the total spin along the 
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direction of the transverse magnetic field b and n+1 = n. 
Our aim is to examine the entanglement between qubits 
at arbitrary sites i,j (i =/= j) in the thermal state 



(2) 



where Z = Tr cxp[— j3H] and T is the temperature (we 
set Boltzmann constant k = 1). This entanglement is de- 
termined by the reduced pair density pij = Tr„ . -{-j-n p(T) 
and can be measured through the concurrence [261 ] 



Ci 



[2X M -tvR}+, R 



1/2- 1/2 
Pij Pij Pij ■ 



(3) 



where [u] + = (u + \u\)/2, Am denotes the largest eigen- 
value of the hermitian matrix R and pij = ^s v i s v ^pijS v i s v ^ is 
the spin flipped density (tri? is the fidelity fil between p^ 
and p^ ) . The entanglement of formation [271 ] of the pair 

is Eij = - J2u=± 1u log 2 q v , where g± = (1± Jl- C?.)/2 



and is just an increasing function of C^- ,with 



Ci- 



1 (0) for a maximally entangled (separable) pair state. 

Since H commutes with S z and is invariant under 
translation and inversion, pij will commute with the pair 
spin component Sfj — sf + s| and its elements will de- 
pend just on the separation \i— j\. Hence, in the standard 
basis of S*j eigenstates, it must be of the form 



Pij 



( vl o \ 

p L a L 

a L p L 

V pi / 



L = 



(4) 



where pi + 2p L + p L = 1, Vl - Vl = 2 ( 



and 



Vl 



(5) 



may rewrite ff exactly as a bilinear form in cj, Cj for 
eac/i value of the spin or fermion number parity 

JV 

P ee exp[iirN] , N = J2 c ] c i = S * + ™/ 2 ■ 



The result for P = er = ±1 i 



is |8J 

i? CT = E & (4 C J ~ I) ~ H 1 _ ^n^l)(ctcj+l + ct +1 Cj) 
3=1 



= ^ ^k{c'lc' k - |), A fe =&-wcosa; fc 



(7) 



where the fermion operators c^ are related to cj by a 
parity dependent discrete Fourier transform 

1 



c 5 = 4t E e -a 



^'4, w k =2-Kk/n, 
{-[in] + Uri, 



(8) 



L ,„ L ,., T [|(n-l)] + i^i} (9) 

with [. . .] denoting integer part. The index k is then 
half-integer (integer) for a = 1 (— 1). 

The 2™ energies are then J^fcgif C^fc — l/2)Afc, where 
iV fc = 0,1 and a = (-1)^*^. Note that the single 
fermion energies A& depend on the global parity a and 
are degenerate (A/- = A_fe) for |fe| ^ 0, n/2. It is also 
apparent from ([7]) that the spectrum of H is independent 
of the sign of 6, and for even n also of the sign of v, as 
coswfc' = — cosu/fc for k' = n/2 — k and kl belongs to the 
same parity as k if n is even. This is also evident from 
([T]), since for even n the sign of v can be inverted by a 
local transformation s*' y — > (— l) 3 s^' y (and that of b by 
sj' z — > — s^' z ). The concurrence (O will then exhibit the 
same properties, depending just on |6| and for even n just 

on \v\ [rn ]. 



Here (O) ee Tr p(T)0 denotes the thermal average of 
O and (sf) = (S z )/n is the intensive magnetization. 
Pij commutes as well with the total spin of the pair 
(S^) 2 — S y • S y , its eigenstates being the standard 
triplet states and singlet Iff), III) and (|f j) ± |IT))/V2, 
with eigenvalues p L , pl ± oil- The pair entanglement 
is obviously driven by the mixing coefficient aL- The 
concurrence ([3]) becomes 



C L = 2 



ax, - 



vLvl 



(0) 



so that pij is entangled if and only if |ax| > \JvlVl- 
This condition also follows from the PPT criterion 12811. 



A. Exact energy levels 

By means of the Jordan- Wigner transformation to 
fermion operators cj = sj" exp[— nrX)fe=i s fe s fc] 0> we 



B. Exact partition function and concurrence 

The partition function Z of the system is to be eval- 
uated in the full grand-canonical (GC) ensemble of the 
fermionic representation. However, due to the parity de- 
pendence of the latter, this requires a number parity pro- 
jected statistics [29]. Z can then be written as a sum of 
partition functions for each parity, 



cr=±l 



CT=±1 



where 4(1 + o~P) is the projector onto parity a and 
Z a v = TrP v e- pH ' = e^ 2 ]J (1 + (-l) v e~ pXk ) , (11) 



fee if „ 



for v = 0, 1. The expectation value of an operator O can 
then be similarly expressed as 



(O) 



\Z-' ]T (Z *(O)S+<7Zf(O)?), (12) 

<7=±1 

(O)^ = (^)- x Tr [pVe-W'O] , u = Q,l. (13) 
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In the case of many-body fcrmion operators, the thermal 
version of Wick's theorem [30| cannot be applied in the 
final average (fl~2|) . but it can be used for evaluating the 
partial averages (JT^) (as p v e~ &H ' = e -PH„+wnN is stm 
the exponential of a one-body operator), in terms of the 
contractions 

g L ee {c\cj)l = ±J2 ^k)Z cos (^) . ( 14 ) 
where {c'\c' k )l = [1 + (-l)^*]" 1 (Eq. Q3J. As s? = 



C Ci 



this leads to 



= So "i (« + 5)(^ + 5))" = 5o 2 -9|- (15) 



Using the identity s+s - = s+[nU + i(4 s fc + s fc 4 )] s 



for i < j, with s+s+ +1 = c]c] +1 , s+s^ = c]c j+ i @|, one 
also obtains 

( S +sp* = iDet(A x ) , (16) 
where Al is the L x L matrix of elements 



(17) 



i.e., Det(Ai) = 2 9l , Det(A 2 ) = 4[ 5l 2 - «? 2 ( ff0 - §)]. All 
terms in (|TJ and ([HJ can then be exactly evaluated. 

In the thermodynamic limit n — ► oo, and for finite 
L -C ?i, we can ignore parity effects and replace sums 
over A: by integrals over w ee cjfe. We can then directly 
employ Wick's theorem in terms of the elements 



9l = (etc 



1 



cos(Lo>) 



-dio . 



(18) 



This leads to (Eq. EJ 

Pt = So " 9l, a L = §Det(A £ ) , (19) 

and = p \ + 1 — 2<7o, where Al is constructed with the 
elements (|18|) . We then obtain the final expression 



|Det(A L )| - 2 v /( 3o 2 -. 9 |)((l-« ?0 )2- ff |) 



(20) 

Note that for T -> 0, Eq. ([11]) yields g L = for 6 > |v| 
and (?l = sin(Lo;)/(L7r) (with go = w / 7r ) for |6| < \v\, 
where cosw = b/\v\. 

When the ground state is non-degenerate, Eqs. (fT9j - 
([20| are also exactly valid for finite n in the T — > limit, 
using the exact contractions 



9i = (4 c j7o = 



o = - cos(Luj k ) , 



(21) 



A' occ. 



where (O)o denotes ground state average and the sum 
runs over the occupied levels (see next section). 



III. RESULTS 
A. Ground state transitions and concurrence 

Let us first describe the behavior in the T — ► limit. 
As [if, N] — 0, the ground state of H can be character- 
ized by the fermion number N, i.e., the total spin com- 
ponent M = N — n/2 in the spin representation. Since 
Xk in ([7j becomes negative for b < v cos Uk , the ground 
state will exhibit n transitions N — > N + 1 as b decreases 
from \v\ to —\v\, starting from N = (the aligned state 
M = -n/2) for 6 > \v\ (X k > V fc) and ending with 
2V = n (M = n/2) for b < -\v\ (X k < V fc). 

For v > 0, the first transition — » 1 occurs at 



i.e., when the lowest negative parity level Ao = b — v 
becomes negative. It represents for b > the entangled- 
separable border at T = 0. For b just below b\ the ground 



state is the one-fermion state 



the W-state (||U ■■■[) + \IU ■ 
exhibits a constant concurrence 



10} = -^£,410), i.e., 
. . i) + . . .)/y/n, which 



Cr 



2/n, (N = l) 



(22) 



for any separation L (Fig. 1). Hence, the transition at 
b = bi is from a fully separable state for b > v (aligned 
state) to a state where any pair is equally entangled. 

Due to the parity dependence of the energy levels, the 
next transition 1 — > 2 (existing for n > 4) does not take 
place when the next A& becomes negative (6 = vcoswjA 
but rather when the lowest a = 1 level crosses the previ- 
ous a = — 1 level, i.e., when 2X±±/2 = Ao, which leads to 
&2 = v{2 cos(7r/n) — 1). In general, for v > the transi- 
tions iV — 1 — > iV occur at 6jy = 2vcoso;fe — 6jv-i) with 
fe = (JV — l)/2, which leads to the critical fields 



b N = v 



cos[{N - \)ir/n 
cos[7r/ (2n)] 



l<N<n, (23) 



i.e., 6at = v(cosuik ~ sinw^ tan(7r/2n). Thus, &at < 6at-i, 
with 6 n _jv+i = — &w and b^ ~ vcosuk for large n. 

Eq. (|22p is valid for 62 < 6 < 61. The exact expression 
for Cl at the other 7V-fermion ground states is given by 
Eq. (|20|) with the elements ([2Tj). which become 



s'm(NLir/n) 
nsin(L7r/n) ' 



(24) 



with g = N/n = \im L ^ 9l- For N = I, g L = 1/nV L 
and Eq. leads to Eq. (|2"2"j) . 

For N > 2, Cl will depend on the separation L, de- 
creasing almost linearly with L for not too small n, as 
seen in Fig. [TJ A series expansion of (|20|) yields the 
initial trend C L w (2iV/n)[l - ixL^J (N 2 - l)/3/n]) for 
A^L -C n. The extent of pairwise entanglement de- 
creases then rapidly as N increases (inset in Fig. [TJ, 
the separation between most distant entangled qubits 
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FIG. 1: (Color online) Top panel: The concurrence Cl as 
a function of site separation L in the 20 different entangled 
ground states of Hamiltonian ^ existing for n = 40 qubits 
and b > 0. Ground states are labelled by the effective fermion 
number N = M + n/2 — 1, . . . , n/2, which can be selected 
by adjusting the magnetic field b (N fermions for &zv+i < b < 
bff)- For N = 1, Cl has the same value for all separations 
(Eq.[22}. The inset depicts the entanglement range L m vs. the 
magnetic field in the same system, which for b > v vanishes at 
T — but remains maximum if < T < T c = T n / 2 (see next 
section). Bottom: Same details for an odd antiferromagnetic 
chain with n = 41 qubits (results for the T — > limit of 
p{T)). For ./V = 1 all pairs are still entangled but Cl decays 
for increasing L (Eq. I25|l . 

being L m « [(n + 1.79)/3.57] for N = 2 and, roughly, 
L m M [(n + 4)/(2iV)] for 2 < N « n/2. 

Just first and second neighbors [L = 1,2) remain en- 
tangled for |6| < 0.65u V n (and |6| < 0.82v (N > n/5) 
for n — > oo) whereas only adjacent pairs [L = 1) re- 
main entangled for |6| < 0.26i; V n ^ 5 (and |6| < 0.5t> 
(iV > for n — ► oo) (for n = 5 second neigh- 

bors are entangled V b > 0). The concurrence of adja- 
cent pairs increases first linearly with N (Ci « 2N/n 
for AT -C n) and becomes maximum for A~ = n/2 
(n > 4), where 51 = l/(nsin(7r/n)) « 1/7T for large n 
and Ci = 2 3 i(l + 5 i) - 1/2 w 0.339. 

For odd n, results for u < must be separately exam- 
ined. The lowest negative parity level is now X±\ n /2] — 
b — \v\ cos(7r/n), so that the first transition occurs at 

= \v\ cos(7t/n) , (v < 0, n odd) , 

with the ground state two-fold degenerate after the tran- 
sition (k = ±[n/2]). The concurrence of the mixture 
\^lk=±[n/2} °f tne tw0 degenerate ground states 



1 



J2 e -^u c t| ) (the T -> limit of p(T)) is 



2cos(L7r/n)/n, (N = 1) 



(25) 



which is again non-zero V L (n is odd) although it now 
decays as L increases (bottom panel in Fig. [T]). For 
L <C n, m 2/n, in agreement with (|22p . whereas for 
most distant qubits (L = [n/2]), = 2 sin[7r/(2n)]/n « 
ir/n 2 . Hence, for large L a significant odd-even difference 
in Cl arises if v < 0, even for large qubit number n, due 
to the ground state degeneracy of the odd system. 

The next transition for v < and n odd occurs when 
A n /2 + K/2-1 = A[ n / 2 ], i.e., at 6^ = |u|[l + cos(27r/n)] - 



b x , and in general at b N — \v \ ^2!k=N/ 
which leads to the smaller critical fields 



V \Y,k=N/2~l C0SUJ k - b N _ v 



b N = 6w cos(7r/n), 1 < N < n , 



(26) 



where 6at are the fields (f2"3"| . Ground states remain two- 
fold degenerate V N ^ 0, n, since there is just one fermion 
in the highest occupied level (k = ±(n — N)/2). 

Eq. ([25]) holds for b^ < b < b± . The expression of C~[ 
for general N in the T — » limit can be similarly obtained 
by using Eqs. PD|) - ([2T1) for each of the degenerate ground 
states and taking then the average. The final result is 



|Re[Det(A-)]| - 2y/(g* - fl £)((l - g f - gj) 



where A L is constructed with the elements 



9l = 9Le 



;L7T /r, 



(27) 



(28) 



with gL given again by (f24|) . For AT = 1 ([27]) leads to 
Eq. The behavior of for A^ > 2 is similar to 

that of Cl (Eq. |2"0"]) , although it is smaller than (due 
to the ground state degeneracy) and its decay with L is 
less linear (see bottom panel). For instance, for L = 1 
Re(Det^4^~) = Det(^4i) cos(7r/n), whence C-f < C±, with 
Cf — > Ci for large n. 

Let us finally mention that for n — ► 00 and nN/n — > w, 
with L finite, Eqs. ([24]) - (128)) coincide both exactly with 
the limit of (TTB]) for T — > 0, where 6at — ► wcosw. 



B. Results for finite temperatures 

Illustrative results for n = 14 — 15 and the thermo- 
dynamic limit n — > 00 are depicted in Figs. [2H2] For T 
close to 0, the concurrence exhibits a stepwise behavior 
in finite chains, in agreement with the T = transi- 
tions previously described, presenting dips at the critical 
fields (f2"3")) - (r2"6"j) due to the ground state degeneracy at 
these points (level crossing) . It is also verified that Cl is 
smaller in odd antiferromagnetic chains, particularly for 
large L close to n/2, in agreement with Eqs. (|2~5|) - (I27|) . 

While at T = there is no entanglement in the ground 
state for b > bi, a fundamental result for T > is that 
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FIG. 2: (Color online) Concurrence (top) and limit temper- 
atures for entanglement Ti,(b) (bottom) for pairs i,i + L as 
a function of the magnetic field, for n = 14 qubits and in 
the thermodynamic limit n — > oo. The concurrence is plotted 
close to the T — > limit (T = 0.005?;). All limit temperatures 
remain constant for b/v — ► oo (see text). Results for even n 
lie mostly above those for n — + oo, particularly for large L 
(where they saturate) and are independent of the sign of v. 



p{T) remains entangled for all fields b > b\ if T is suffi- 
ciently low, leading to a small but non-zero concurrence 
Cl for any separation L if < T < Tl(o). Moreover, 
the limit temperature ?l(6) approaches a non-zero limit 
Tl for b — ► oo V L, being practically constant for b > \v\ 
(and V6 if L = 1). This behavior applies for any n, in- 
cluding n — > oo as well as the special case v < and n 
odd, as seen in the lower panels of Figs. [2][3l 

In order to rigorously prove the previous behavior, we 
note that for b — \v\ > kT (e _/3 ( 6 ~H) ^ i) : we may keep 
just zero, one and two fermion states in exp[— f3H], i.e., 



k<k'eK + 



-/3(A fc +V)l 



and similarly for p{T). This leads to «i « e (f3v) 

and pt ~ e _2/3& [/(^ 2 (/3w) — I^ 2 (f3v)] up to lowest order 
in ', where 



Jt(0«) = - E e^ cos ^cos(Lo; fe ) 



k£K ± 



Hence, up to first order in e wc obtain 



CY « 2e 



-/3b 



(29) 



(30) 



Thus, as b increases the concurrence decreases exponen- 
tially with the field when it is positive, but the limit tem- 
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FIG. 3: (Color online) Same details as Fig.[2]for n = 15 qubits 
and v < 0. Results now lie mostly below those for n — *• oo 
(which are the same as in Fig. [2} and do not saturate for large 
L. Results for n = 15 and v > are similar to those of fig. [2] 

perature T^(b) becomes constant, as the entanglement 
condition Cl > becomes 6-independent, i.e., 



-((3v) > Jl+ 2 (0v)-lf(l3v) 



(31) 



Eq. (|3ip is always satisfied for sufficiently small but pos- 
itive T, for any distance L, ensuring a non-zero con- 
currence and limit temperature for any b > \v\. 
This is easy to prove for v > 0, where for T — > + , 
« e^/n > I + (/3w) « 2e /3 " cos ( 7r / Tl Vn. It also 
holds for v < (n odd), since in this case, for T — > + , 
!£" 03u) « e /5WI^(^/™) cos(Ltt /n)/n, whereas the r.h.s of 
plj) becomes » ™*< 2 (fA0 sin(Z,7r/ra)/ra < (/3t>). 

In the thermodynamic limit n — > oo, and for finite 
L <C n, we may neglect parity effects and just replace 



Il(f3v) [ e^ VCOSLJ cos(Lu>)du = I L {(3v) 

n Jo 



(32) 



where Il(x) is the modified Bessel function of the first 
kind [Il{x) « e x [l + (1 - 4L 2 )/8a;]/\/2^ : for a; -> oo, 
with I L (-x) = (-) L / L (a;))- Eq. ^3DJ becomes then 
identical with the result obtained from Eqs. HHJ-ffJOl) 
(9L e-P b I L (/3v) for b-\v\ » T, with Det(A L ) -» 2 ffL ). 
Eq. (j3"Tj) becomes then 



V2Jl(/?M) > /o(/3«) 



(33) 



which is again always satisfied for sufficiently low T V L. 
The limit temperatures Tl = Tl{oq) are then determined 
for n — ► oo by the equation V^l^M) = -^o(/3|w|), which 
leads to T x w 0.486|v|, T 2 w 0.16|d and 



T L « |u|ln2/L 2 



(34) 
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FIG. 4: (Color online) Top: The limit temperatures Tl for 
entanglement between pairs for large magnetic fields b ^> v 
as a function of separation L, for n — 14 and n = 40 qubits, 
determined by Eq. (I31|l . Inset: The ratio between Tl and the 
value in the thermodynamic limit T£°, determined by Eqs. 
(|33l) -(l34l. Tl deviates from T£° for L > n/4, approaching, 
for L — > n/2, Eq. (|36[l for n even or n odd and v > 0, and Eq. 
(|37l) for odd and v < 0. Bottom: The limit temperature for 
the most distant qubits (L = [n/2]), showing the odd-even 
staggering arising for v < (dashed line). The upper (lower) 
dotted line depicts the result of Eq. J36J ((EH). 



for large L (as I L {x)/I {x) w e - i2 /(2*) for x > L 2 ). 
Thus, T L (b) decreases as the inverse square of the 
pair distance L for large b. The maximum value at- 
tained by Cl for b > \v\ becomes nevertheless small 
and decays exponentially with both b and L 2 (Cl ~ 
e -(b/\v\-l)L*/tf( t y L for T = \ v \ t jL 2 < T L , with f(t) = 
v /2i/7r[e"'/ 2 -Vl - e"']). Eq. ^ also indicates roughly 
the value of Tl(6) at the critical region b « \v\, since it 
is almost constant for b > |w| (Figs. HE]). 

On the other hand, for large L « n/2 the projected 
expression (|30|) is required even /or /ar^e n. For instance, 
for even n and L = n/2, cos(Lujk) = ((— 1) ) for k half- 
integer (integer). Hence, in this case I^/2^P V ) ~ ^' whil 6 
for v > and large n, 



6*2(4) (e 



-2/3vir 2 /n 2 



)/n, (35) 



1-^/2 (e 2 ( U )^2Er=i 



/2 



after replacing coswfe 

6» 4 (u) = 1 + 2^]^ 1 (-l) fe u feZ , denote the Elliptic Thcta 
functions). These results also hold approximately for 
large odd n and L = [n/2]. Eq. (j3"Tj) becomes then 
I~, 2 (f3v) > Iq((3v), and since 9 2 (u) — 6* 4 (u) for u = e _7r , 
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FIG. 5: (Color online) The tr/ermal behavior of Cl for n = 
40 qubits near the transition at b = v. Limit temperatures 
remain stable at the transition, indicating the reentry of Cl 
for T > V L for b > v. The inset is an enlargement of the low 
T region, showing the accumulation of the limit temperatures 
for large L at a value close to that given by Eq. (13611 . 



it leads to the limit temperature 

T[ n / 2 ] ~ 2tto/ n 2 . 



(36) 



for the most distant pairs and v > 0. It is greater than 
Eq. ([33]) for L = n/2 by a factor vr/(21n2) ps 2.27. 

Eq. (|3l)l) does not hold for w < if n is odd. In this 
case we may directly employ the asymptotic expression 
of Eq. (|31[) for T — > + , which for large n yields 



n/2] 



\v\ir 2 



In'- ln[2\/2n/ 7T 



(v < 0, n odd) . (37) 



Thus, in this case there is an additional logarithmic factor 
in the denominator, which makes T , [„y 2 ] lower than Eq. 
([55)1 and also Eq. for L — n/2, originating an odd- 
even staggering of T„/ 2 if v < 0. 

The behavior of Tl is depicted in Fig. [4j It is seen that 
for L > n/4, it deviates from the 1/L 2 law given by Eq. 
([3~4")l . approaching the values given by Eqs. ([55)) or ([3T)) 
for L n/2. Fig. [5] depicts the typical thermal behavior 
of Cl for u > near the transition at b = b%. For b <b\ 
there is entanglement between all pairs if T is lower than 
a certain temperature, given approximately by Eq. (|36[) . 
It also shows the reentry of Cl for T > for b > v, which 
is quite prominent for low L. 

Finally Fig. [6J depicts the typical behavior with the 
qubit number n of the concurrence and limit tempera- 
ture. We have chosen a separation L = 2. Although 
the thermodynamic limit is on the average rapidly ap- 
proached, the stepwise behavior of Cl at low T « Q.Olv 
remains visible even for n = 40, and deviations in the 
limit temperature can be significant at the onset. They 
are as well significant for small n < 10. 

An interesting feature is that the slope of Tl(&) can 
be negative in this region, a fact already seen in Fig. [2] 
for n = 14, and visible here for n = 8 and n = 20. This 
occurs when the value of the onset field b c for finite n 
(which for L — 2 corresponds to bi, 63, 67 and &14 for n — 
6, 8, 20 and 40) lies above the value for n — > 00, as occurs 
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FIG. 6: (Color online) Concurrence (top) and limit temper- 
atures for entanglement (bottom) for separation L — 2 and 
different n, as a function of the magnetic field. The con- 
currence is plotted at T = O.Olv. Dotted lines depict the 
thermodynamic limit. 

for n — 8, 20. In these cases there is a small field interval 
below b c where entanglement between second neighbors 
exists only above a threshold temperature T L (b) > 0, up 
to the higher limit temperature T^ib). 

A final comment is that we have checked all expressions 
by comparison with calculations for low n < 10 based on 
the direct diagonalization of H. In particular, for n = 2, 
the entanglement condition (pTTj) becomes exact V b, as 



in this case there are just one and two-fermion states, 
reducing to sinh(/3w) > 1 and leading to the known limit 
temperature Ti = v/ ln(l + \/2) 0- 



IV. CONCLUSIONS 

We have provided an exact analytic treatment of 
the entanglement between arbitrary pairs in cyclic XX 
chains in the presence of a transverse magnetic field, valid 
both at zero and finite temperatures and for any qubit 
number n. We have shown that in spite of its simplicity, 
this system exhibits very interesting features such as a 
discrete set of [n/2] different entangled ground states at 
T = (and b > 0), which can be easily selected by ad- 
justing the magnetic field across the critical values (|23|) or 
(|26| . and which develop increasing entanglement ranges, 
reaching always full range (all pairs entangled) in an in- 
terval &2 < b < b±, even for odd antiferromagnetic chains. 

Moreover, while at T = the ground state is fully 
separable for b > b\, we have rigorously proved that for 
T > there is a small but non-zero entanglement be- 
tween any pair for all fields b > bi if T is sufficiently low, 
which decays exponentially with the field and with the 
square of the separation L. Limit temperatures Tl are 
roughly independent of b for b > v and decay as L~ 2 for 
L < n/4, but tend to saturate at 2~[ n /2] (Eq-EH]) for most 
distant pairs (L ps n/2) if n is even or v > 0. We have 
also shown that due to degeneracy of the ground state, 
pairwise entanglement in odd antiferromagnetic chains 
is weaker, particularly for distant pairs, where odd-even 
effects in the concurrence and Tl subsist for all n. 

NC and RR acknowledge support of CONICET and 
CIC of Argentina. 
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